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ABSTRACT 

We address the issue of accuracy in recovering density profiles from observations of rotation 
curves of galaxies. We "observe" and analyze our models in much the same way as observers do 
the real galaxies. Our models include stellar disks, disks with bars, and small bulges. We find 
that the tilted ring model analysis produces an underestimate of the central rotational velocity. 
In some cases the galaxy halo density profile seems to have a flat core, while in reality it does not. 
We identify three effects, which explain the systematic biases: (1) inclination (2), small bulge, 
and (3) bar. Inclination effects are due to finite thickness of disk, bar, or bulge. Admixture of 
a non-rotating bulge component reduces the rotation velocity. A small (200-500 pc) bulge may 
be overlooked leading to systematic bias even on relatively large ~ 1 kpc distances. In the case 
of a disk with a bar, the underestimate of the circular velocity is larger due to a combination 
of non-circular motions and random velocities. The effect of the bar depends on the angle that 
the bar makes with the line of sight. Signatures of bars can be difficult to detect in the surface 
brightness profiles of the model galaxies. The variation of inclination angle and isophote position 
angle with radius are more reliable indicators of bar presence than the surface brightness profiles. 
The systematic biases in the central ~ Ikpc of galaxies are not large. Each effect separately gives 
typically a few km/serror, but the effects add up. In some cases the error in circular velocity 
was a factor of two, but typically we get about 20 percent effect. The result is the false inference 
that the density profile of the halo flattens in the central parts. Our observations of real galaxies 
show that for a large fraction of galaxies the velocity of gas rotation (as measured by emission 
lines) is very close to the rotation of stellar component (as measured by absorption lines). This 
implies that the systematic effects discussed in this paper are also applicable both for the stars 
and emission-line gas. 

Subject headings: galaxies: halos — galaxies: local group — galaxies: dwarf — galaxies: kinematics and 
dynamics — cosmology: dark matter — 
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1. Introduction 

The cosmological cold dark matter (CDM) 
model and its present variant with cosmologi- 
cal constant (ACDM) have been very successful 
in accounting for the properties of galaxies and 
their distribution in space. A current problem for 
CDM lies in observations of dwarf and low sur- 
face brightness galaxies. The rotation curves and 
derived density profiles of these galaxies appear 
to be in conflict with the CDM model (Moore 
1994; Floras & Primack 1994; Burkert 1995; 
Hernandez & Gilmore 1998). The problem arises 
because numerical simulations of the formation 
of halos in the CDM model reveal steeply cusped 
density profiles while observations favor shallow 
(core) profiles. The theory predicts that the mass 
density in the inner parts of the dark matter ha- 
los is best described as an r~", a = 1 power law 
(Navarro, Frenk, & White 1997; Power et al. 
2003; Colin, Klypin, Valenzuela, Gottkx;ber 2004) 
or possibly by a power-law with a slightly steeper 
a =1.2- 1.5 slope (Moore et al. 1999; Ghigna et 
al. 2000; Klypin, Kravtsov, Bullock, & Primack 
2001; Reed et al. 1999; Tasitsiomi, Kravtsov, Got- 
tloeber, & Klypin 2004). The observed rotation 
curves of a significant fraction of dwarf galaxies, 
however, suggest that their cores do not have steep 
density cusps (de Blok & McGaugh 1997; Blais- 
OucUcttc, Amram & Carignan 2001; dc Blok, 
McGaugh, Bosma, & Rubin 2001; McGaugh, Ru- 
bin, fc de Blok 2001; Bolatto, Simon, Leroy, & 
Blitz 2002; Weldrake, dc Blok & Walter 2002; 
Simon, Bolatto, Leroy, & Blitz 2002). Note that 
this is not true for some galaxies and that some 
observational results are still consistent with cos- 
mological predictions (van den Bosch et al. 1999; 
Swaters, Madore, & Trewhella 2000; van den 
Bosch & Swaters 2001). McGaugh, Rubin, & de 
Blok (2001) have analyzed Hq observations of 
50 low surface brightness galaxies of which 36 are 
of sufficient quality that rotation curves can be 
derived. The spatial resolution is improved by an 
order of magnitude over old 21 cm data. In a 
companion paper de Blok, McGaugh, Bosma, & 
Rubin (2001) conclude that at small (0.5-1 kpc) 
radii, the mass density distribution is dominated 
by a nearly constant density core. The distri- 
bution of inner slopes is strongly peaked around 
a = —0.2 with clear contradiction with the CDM 
predictions. Thus, there appears a controversy 



concerning the density distribution of matter in 
the inner parts of dwarf galaxies. 

We propose a somewhat different approach to 
this problem which we believe can help to clarify 
these debates. We construct realistic models of 
galaxies and "observe" them as real galaxies using 
the same techniques or algorithms that observers 
apply to real data. We study the details of the 
analysis of observations and we find discrepancies 
between recovered and true properties of the sim- 
ulated galaxies. We investigate projection effects 
and the effects of bars and bulges on the observed 
rotation curves. We show that these effects are 
often far from negligible. 

Some of the effects which we study, have al- 
ready been discussed. Athanassoula (1984) dis- 
cussed the biases in the observed rotation curves of 
high surface brightness galaxies created by strong 
bars. More recently, Koda & Wada (2002) have 
investigated the effects of a weak bar on mass de- 
termination in the central parts of spiral galaxies. 
They examine rotation curves of edge-on model 
galaxies observed from different viewing angles. 
The models had a two-dimensional solid rotating 
barred potential. Motion of gas in this potential 
was followed using an SPH code. Koda & Wada 
(2002) find that the mass estimated from rotation 
curves is larger than the true mass by a factor 
of 5 for 15% of the viewing angles. Overestima- 
tion in mass occurs more frequently than under- 
estimation. Koda & Wada (2002) focus on rota- 
tion curves defined as the highest velocity envelope 
of position-velocity diagrams. When the rotation 
curve is defined as the intensity weighted mean of 
position- velocity diagrams, the rotation velocity is 
underestimated for most of the view angles. Koda 
& Wada (2002) conclusions are only for edge-on 
galaxies; the situation with less inclined galaxies 
could be more complex. Another important aspect 
studied by Wada, Meurer & Norman (2002) and 
Wada & Koda (2001) is the gravity driven tur- 
bulence induced by a weak bar in a clumpy and 
filamentary gaseous media. If was found that the 
velocity dispersion of a cold neutral gas as well 
as of an ionized component is significantly lager 
than naively expected if the distribution is clumpy 
(Wada, Meurer & Norman 2002) or filamentary 
(Wada & Koda 2001). These results challenge 
the common expectation that the gas component 
in a galaxy is a cold medium with low velocity 
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dispersion. 

The effect of bulges on the observed kinemat- 
ics of galaxies has its own history (e.g., Fillmore, 
Boroson, & Dressier 1986; Kormendy & Westpfahl 
1989; Bertola et al. 1995; Vega Beltran et al. 2001; 
Funes et al. 2002; Pizzella, Corsini, Vcga-Beltran 
& Bertola 2004). An interesting issnc here is the 
difference between the rotation of stars and rota- 
tion of gas in the central 1-2 kpc region. Because 
our simulations have only the stellar component 
(no gas), it is important to know how different are 
rotation of gas and of stars. Naively, one expects 
that g "colder" component should rotate 

faster. The reality is more complicated. A de- 
tailed discussion of this issue is given in §2. We 
also use our own observations to measure rotation 
velocities in stellar and gas components. To sum- 
marize, like other workers in this field, we find that 
often in the central region the rotation of gas (as 
measured by Ha and Nil emission) is not very 
different from the rotation of the stellar compo- 
nent. 

Some other systematic effects were considered 
by Swaters, Madore, van den Bosch, & Balcells 
(2003) who discuss effects of seeing, finite slit 
width, and slit misalignment for ideal thin expo- 
nential disks. It is found that all effects systemat- 
ically underestimate the central slope. The mag- 
nitude of the bias depends on the distance to the 
galaxy and on slit misalignment, de Blok, Bosma, 
& McGaugh (2003) considered the effects of an 
offset between optical and kincmatical centers, as 
well as a very simplistic attempt to quantify the ef- 
fects of non-circular motions. Both effects tend to 
lower the estimated density profile slope. Hayashi 
et al. (2003) discuss effects of triaxial halos on the 
rotation of gas in CDM simulations, which lead to 
non-circular motions. 

This paper is organized as follows. In §2 we 
present results of our observations of rotation ve- 
locities as measured by gas and stars. In §3 we 
discuss different observations, focusing on system- 
atic biases, which may affect rotation curves. We 
also present resiilts of our observations of rotation 
of gas and stars. In §4 we present and describe 
our models. Simulations are discussed in §5. We 
have constructed four detailed galaxy models and 
run them as N-body simulations. In §6 we present 
our methods of analysis of the models and in §7 we 
show our results. In §8 we present our conclusions. 



2. Observational measurements of rota- 
tion curves: gas vs stars 

Most observational results which address the is- 
sue of cusps and cores, provide rotation for gas 
as measured by either Ha emission or by 21 cm 
line. The rotation of the stellar component is 
not studied because it is difficult to measure for 
dwarf and LSB galaxies. Nevertheless, there are 
observational results for stellar rotation curves for 
a substantial number of galaxies (e.g., Fillmore, 
Boroson, & Dressier 1986; Kormendy & Westp- 
fahl 1989; Bertola et al. 1995; Vega Beltran et al. 
2001; Funes et al. 2002; Pizzella, Corsini, Vega- 
Beltran & Bertola 2004). The observed set of 
galaxies cover a wide range of morphological types 
(from SO's to Sc's) and magnitudes. Most of the 
galaxies are normal high surface brightness galax- 
ies, but some are similar to dwarfs studied for the 
"cusps and cores" problem. 

Let's start with what one naively should ex- 
pect. Close to the center stars are expected to 
have relatively large random (non-circular) veloc- 
ities. For example, if the stars are in a disk with 
a constant height, the velocities perpendicular to 
the plane of the disk should increase exponentially 
as one gets close to the center in order to react to 
the exponentially increasing disk density. Other 
random velocity components are coupled with the 
vertical velocities with the radial rms velocity be- 
ing typically twice larger than the vertical velocity. 
Non-circular stellar velocities may be even larger 
if there is a bulge. Large non-circular velocities re- 
sult in significant asymmetric drift. In turn, this 
would imply that rotational stellar velocity falls 
significantly below the circular velocity. In the 
case of gas, one naively expects that rotation is 
faster because gas cools and settles in the plane of 
the disk. 

The problem with this picture, which is shared 
by many astronomers, is that it is too naive. It 
has some justification in the outer parts of disk 
galaxies. Indeed, in our own galaxy at the solar 
neighborhood molecular gas has a smaller height 
than stars. Thus, gas has smaller non-circular ve- 
locities and must rotate faster. Yoimg stars, which 
have formed relatively recently from the gas, have 
smaller asymmetric drift than the older stars. This 
picture of a cold thin gaseous disk inside a thicker 
slower rotating disk or a bulge is naively applied 
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to the central parts galaxies. Unfortunately, close 
to the galactic center the naive picture has lit- 
tle justification. For example, the disk cannot be 
considered as thin because at distances 0.3-0.5 kpc 
from the center the disk thickness is comparable 
with the distance to the center. Gas motions also 
increase substantially (Vega Bcltran ct al. 2001; 
Funes et al. 2002; Pizzella, Corsini, Vega-Beltran 
& Bertola 2004) 

Rotation curves of ionized gas (Ha) as well 
as rotation and velocity dispersion curves for the 
stellar component (absorption lines) were studied 
by Fillmore, Boroson, & Dressier (1986). It was 
found that for some galaxies the stellar component 
rotates faster than the ionized gas in the central 
region. For other galaxies the rotation of both 
components is comparable. In particular, for the 
Sab galaxy NGC 4569, it was found that at 300 pc 
the ionized gas rotates with a speed of ~50 km/s, 
while the rotation detected in the stellar compo- 
nent is ~100 km/s and the stellar velocity disper- 
sion is ~150 km/s. The situation is similar inside 
500 pc. This led Fillmore, Boroson, & Dressier 
(1986) to suggest that gas rotation might not rep- 
resent the true circular velocity. NGC 5879, an 
Sb galaxy, shows the same rotation speed for ion- 
ized gas and for the stellar component up to the 
distance of 1 kpc. However, in addition to the 
rotation velocity the stellar component has a ran- 
dom component. At 200 pc both stars and ionized 
gas have a rotation speed of 40 km/s, but the stel- 
lar component also shows a velocity dispersion of 
70 km/ s. Based on these observations as well as on 
their mass models, Fillmore, Boroson, & Dressier 
(1986) suggest that ionized gas in these galaxies 
is likely to be supported not only by rotation but 
also by pressure. 

For NGC 4594 (M104) Kormcndy & Wcstp- 
fahl (1989) measured the rotation of ionized gas 
(HII regions) and inter-arm gas using Ha obser- 
vations. The rotation of the stellar component 
was measured using absorption lines. It was found 
that in the central region ionized gas rotates with 
a 200 km/s, while the stellar component rotates 
faster: its rotation velocity is about 300 km/s. 
Kormendy & Wcstpfahl (1989) suggest that ion- 
ized gas in M104 has smaller rotation because it 
is either shed by the bulge stars or is somehow 
related with the x-ray halo. The situation for 
galaxies with small bulges is not clear. Gas in 



small (few hundred parsec) bulges with an effec- 
tive radius comparable to the disk characteristic 
height is probably dense enough to harbor star 
formation. A coupling between young stars and 
ionized gas could produce non-circular motions in 
the gas kinematics and lead to an underestima- 
tion of the central mass content. Other possi- 
bilities are discussed and explored by Cinzano et 
al. (1999). Among these, is a model where gas is 
made of clouds. A similar situation was discussed 
also by Wada & Koda (2001) and Wada, Meurer 
& Norman (2002). The clumpy and filamentary 
structure of a gaseous component considerably in- 
creases the asymmetric drift too. Wada, Meurer 
& Norman (2002) discuss the case of NGC2915, a 
blue compact dwarf galaxy, that shows a bar and 
a large velocity dispersion in the H-I kinematics. 

There are no strict rules as far as rotation of 
gas and stars in the central parts of galaxies are 
concerned. In some cases the gas rotates faster 
as naively expected (e.g., Sb galaxies NGC772 
and NGC 7782 (PignateUi et al. 2001)). In many 
cases, the rotation of gas and stars in the central 
1-2 kpc is nearly the same (e.g., Sc NGC 5530 
and Sb NGC 615 (Pizzella, Corsini, Vega-Beltran 
& Bertola 2004)). Nevertheless, there is a gen- 
eral trend that late spirals with small bulges have 
gaseous and stellar components rotating with al- 
most the same velocities (e.g., Cinzano et al. 1999; 
Vega Beltran et al. 2001). 

The observational data used in our paper are 
from a program aimed at understanding line 
strengths and line strength gradients in bulges. 
Long-slit spectra were obtained along the ma- 
jor axes of 37 spiral galaxies using the Double 
Imaging Spectrograph (DIS) on the ARC 3.5m 
telescope at Apache Point Observatory. The spec- 
tra cover a large wavelength range (4000 to 7500 
A or larger) at the cost of spectral resolution (6-8 
A per pixel). The galaxies span a wide range of 
inclination angles and Hubble types. Some of the 
inclinations are too low to measure rotation and 
the condition that Ha is present in emission elim- 
inates several, mostly early-type, galaxies. From 
the sample galaxies with different morphological 
types and inclinations we selected all 8 galaxies, 
that had inclinations in the range 45-80 degrees 
and had Ha in emission at the very center and 
through the most of galaxy. We did not specif- 
ically target late type galaxies, but, as one may 
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expect, the condition that Ha is in emission at 
the galactic center ehminates most early types. 

Spectra extracted at different positions along 
the sht were cross-correlated with the spectrum at 
the galaxy center. For stellar rotation, the cross- 
correlation was performed within the wavelength 
range 4100 to 5400 A after masking out emis- 
sion features and sky lines. For 6 galaxies, rota- 
tion in the gas component was measured by cross 
correlating within the wavelength ranges 6550 to 
6600 A to include Ha and Nil A6583. For the 
remaining two galaxies (NGC 6368 and IC 1029) 
we used SIX doublet (6548-6583 A ) because it pro- 
vided higher signal-to- noise ratio than the Ha. A 
detailed description of the observations and data 
analysis will be presented in an upcoming paper 
(Moorthy & Holtzman, in preparation). Table 1 
gives properties of these galaxies. 

Rotation curves of the galaxies are shown in 

Figure 1. Surface brightness profiles along the 
spectrograph slit are shown in Figure 2. Both the 
velocities and the surface brightnesses are mea- 
sured along the major axis of the galaxies. For 
barred galaxies they are long the bars. A large 
fraction (but not all) galaxies in our sample show 
stellar rotation close to that of the gas. All three 
barred galaxy in our sample, have complicated ve- 
locity patterns. In case of NGC5719 there is lit- 
tle stellar rotation along the bar. Yet, the bar in 
the galaxy is quite extreme: more than 10 kpc. 
This galaxy also presents a peculiar morphology 
with a single strong dust lane oriented in different 
direction with respect to the galaxy major axis. 
A large bulge dominates the central light distri- 
bution. Based on spectroscopic observations in 
visible, near and far infrared light as well as in 
radio wavelengths Hill, et al. (1999) and Hill 
et al. (2001) conclude that NGC5719 is a com- 
posite of an AGN and a starburst galaxy. This 
evidence combined with the large size of the bar 
and the existence of a companion NGC 5713, sug- 
gests that NGC5719 is suffering a strong interac- 
tion. Hunter et al. (2002) recently presented the 
case of NGC4449. This galaxy is also strongly 
barred and shows no rotation in the stellar com- 
ponent. Hunter et al. (2002) propose a model 
where the orientation of a gaseous younger disk 
and the stellar component are different as a re- 
sult of an interaction. The model introduced by 
Hunter et al. (2002) supports our interpretation 



of NGC5719, however a detailed modeling of this 
peculiar galaxy is beyond the scope of this paper. 
The other barred galaxies in our sample show large 
and strong bars too. 

We clearly find a tendency that stars in galaxies 
with smaller bulges rotate as fast as the gas. This 
is compatible with previous results (e.g., Cinzano 
et al. 1999; Vega Beltran et al. 2001). There are 
important conclusions, which we get form these 
results. (1) Because the stars and the gas move 
in the same gravitational potential, observed Ha 
rotation curves must experience the same system- 
atic biases as the stellar rotation curves. (2) We 
can use stellar-dynamic simulations to study those 
effects. (3) The environment of a galaxy should 
be taken in account in the interpretation of the 
kinematics, specifically for the case of strongly in- 
teracting galaxies. 

3. Observational measurements of rota- 
tion curves: systematic effects 

For most galaxies corrections for finite disk 
thickness seems to be quite small. Indeed, at few 
kiloparscc distance from the galactic center correc- 
tions for disk height of few hundred parsec are not 
significant. As far as we know, no correction for fi- 
nite disk thickness was ever made. Is this justified 
for central regions (0.5-lkpc) of dwarf galaxies? It 
is far from obvious. Observations of edge-on galax- 
ies indicate that as the function of decreasing cir- 
cular velocity the disk height decreases slower than 
the exponential disk scale length (Kregel, van der 
Kruit, & de Grijs 2002; Dalcanton 2003). In other 
words, dwarf galaxies are relatively thicker than 
giants. Using I-band photometry for 34 edge-on 
spiral galaxies Kregel, van der Kruit, & de Grijs 
(2002) give estimates of the exponential heights. 
For 6 galaxies with maximum circular velocity be- 
low 100 km/s their results give the average height 
of 440 pc. One of the thinnest galaxies is UGC 
7321 (Matthews 2000) with an exponential height 
of 150 pc, which is about what we use for our 
model IV "Disk" . Interestingly enough Matthews 
& Wood (2003) find that the neutral hydrogen has 
about the same height as the stellar component. 
If the galaxies, which are used to study rotation 
curves and the dark matter density have the same 
heights as those observed edge-on, it is clear that 
in the central sub-kiloparsec region, the disks can- 
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not be treated as thin at all. 

The effect of finite disk thickness increases with 
increasing galaxy inclination. As part of our study 
we have examined the data of McGaugh, Rubin, & 
de Blok (2001) and de Blok, McGaugh, Bosma, 
& Rubin (2001) to see if we can find evidence 
of observational effects involving inclination and 
galaxy morphology. McGaugh, Rubin, & dc Blok 

(2001) list the inclinations for a sample of 50 galax- 
ies which they observed in Ha. de Blok, McGaugh, 
Bosma, & Rubin (2001) list values for the inner 
power law slope (a) of p(r) for 48 galaxies. There 
is a partial overlap with the sample of de Blok, 
McGaugh, & Rubin (2001), other galaxies being 
culled from data published by de Blok, & Bosma 

(2002) , Swaters, Madore, & Trewhella (2000) and 
Verheijen (1997). We have combined these data 
to produce a list of 43 galaxies. These data are 
shown in Figure 3. If the estimates of the inner 
density slope a had no biases, there would be no 
dependence of a with the inclination. This clearly 
is not so: the top left corner of the plot does not 
have points. 

Also note that galaxies identified as barred 
galaxies tend to have larger values of a than the 
average. All but one have a > —0.25. We do not 
have morphological information for all galaxies in 
the sample. The data in Figure 3 suggest that 
galaxies with high inclinations have large a val- 
ues. Galaxies with inclinations > 60° have a val- 
ues greater than -1. The evidence suggests that 
projection effects and the presence of bars influ- 
ences must be taken into account when determin- 
ing density distributions from rotation curves. 

Recent observations of galaxies suggest that 
galactic bars occur more frequently than had been 
previously suspected. Eskridge et al. (2002) have 
presented BVRJHK imaging of a sample of 205 
bright nearby spiral galaxies. An examination of 
Table I in their paper reveals that in the B-band 
103 of the galaxies arc classified as barred spirals 
whereas in the H band 140 or ~70% of the sample 
are classified as barred spirals. The conclusions 
we draw concerning the effects of bars on the in- 
terpretation of spiral galaxy rotation curves may 
well be relevant to the majority of spiral galaxies. 

Dwarf and LSB galaxies used for analysis of ro- 
tation curves and the dark matter content typi- 
cally are selected to be "bulge-less" : types Sc-Sd. 
Yet, it appears that many of them still have small 



bulges. Boker, Stanek & van der Marel (2003) 
have obtained HST WFPC2 images of a sample of 
19 nearby spiral galaxies of late Hubble type. For 
these galaxies bulges arc not evident in the ground 
based optical data for these galaxies. The HST 
images were combined with the existing ground 
based optical images to obtain surface brightness 
profiles. The conventional wisdom has it that the 
sequence from early to late-type galaxies ranges 
from bulge dominated to disk dominated galax- 
ies. The reality seems to be that only ~30% of 
the galaxies in the sample are pure disk galax- 
ies (Boker, Stanek & van der Marel 2003). The 
remainder are well fitted by adding a bulge com- 
ponent. Some of these bulges are small, extend- 
ing out to ^ 300 pc. Larger bulges extend out to 
2 kpc. As in the case of bars, the evidence suggests 
that the presence of bulges in spiral galaxies must 
be taken into account in the investigation of rota- 
tion curves of spiral galaxies. Galaz, Dalcanton, 
Infante, & Treister (2002) confirm these obser- 
vations in a near infrared study of 88 low surface 
brightness galaxies and normal spirals. They note 
that many of the low surface brightness galaxies 
classified as late types based on optical B-band 
images show a dramatic bulge in the near-IR. The 
data support the conclusion of Beijersbergen, de 
Blok & van der Hulst (1999) that many low sur- 
face brightness galaxies could be regarded as high 
surface brightness bulges embedded in low surface 
brightness disks. Since the surface brightness in 
the near-IR closely reflects the stellar surface den- 
sity it is important to take bulges into account 
when modeling the rotation curves of low surface 
brightness galaxies. 

Note that the bulges are really small. A 300 pc 
bulge is hardly thicker than the disk. Physical 
conditions in such small bulges are not much dif- 
ferent from those in the disk. Stars can form in 
those bulges. 

The next important point to keep in mind is 
that of accuracy. How accurately can one de- 
termine the inclination of a galaxy? We refer 
to NCG 3109, a Magellanic type spiral that we 
have attempted to model in some detail (Valen- 
zuela, Rhee & Klypin 2004). Jobin & Carig- 
nan (1990) quote several values of inclination for 
this galaxy of 80° ± 2° based on B-band photom- 
etry. The same authors flnd a value of 75° ± 2° 
based on photometry of an I-band plate taken us- 
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ing the U.K. Schmidt telescope. In a more recent 
study, Blais-Ouellette, Amram & Carignan (2001) 
present CFHT Fabry-Perot H-a observations of 
NCG 3109. These data show an inclination be- 
tween 85° and 90° with one inner point at 74°. 
Although the formal uncertainty in inclination in 
a given waveband may be as small as 2°, a more 
realistic value is at least 5° for the uncertainty in 
a galaxy's inclination. Barnes & Sellwood (2003)) 
support this conclusion. Based on a study of 74 
spiral galaxies they find that systematic uncertain- 
ties in position angle and inclination due to the 
presence of nonaxisymmetric structure are of or- 
der ~ 5°. These estimates are of course lower lim- 
its on the real errors. Weldrake, de Blok & Walter 
(2002) publish inclinations and position angles of 
the local group dwarf NGC 6822 as a function of 
radius based on radio data obtained with the Aus- 
tralia Telescope Compact Array. The values for 
the inclination angle of this galaxy as a function 
of radius fluctuate between ~ 80° and 20°! In their 
analysis, a constant value 60° is adopted without 
further explanation. The position angle for this 
galaxy varies between 100° and 140° according to 
these data. 

To summarize, there are three key points that 
have guided our analysis. Firstly, the central den- 
sity profile obtained for a given galaxy appears to 
be influenced by the inclination of the galaxy, for 
inclinations greater than 60° at least (it is difiicult 
to determine the inclination of these dwarf galax- 
ies to better than several degrees; ~ 5°). Secondly, 
the data seem to indicate that galaxies with bars 
have densities that flatten in their centers when 
the standard analysis is performed. Errors of sev- 
eral degrees (~ 5°) can occur when determining 
inclinations. Thirdly, given increasing evidence for 
bulges in spiral galaxies of all morphological types, 
it is important to examine the biases that bulges 
(as well as bars) introduce into the determination 
of central density profiles. 

4. Models 

We use A^-body simulations to produce mod- 
els of galaxies. The models include realistic dark 
matter and stellar component, but do not have 
gas. Gas would significantly complicate numer- 
ical models making predictions less certain. In 
dwarf galaxies gas is very patchy and often has 



significant non-circular motions (e.g., Bolatto, Si- 
mon, Leroy, & Blitz 2002; Simon, Bolatto, Leroy, 
& Blitz 2002). The stellar component has its own 
complications as we find below, but those are rel- 
atively simple and easy to interpret as compared 
with the gas motions. In dwarf galaxies, the gas 
has all the complications of stars and in addition 
it brings its own problems. Studying the stellar 
component is definitely a good starting point. 

We study the kinematics and spatial distribu- 
tion of stellar component in four models. Parame- 
ters of the models are given in Table 2. All the 
models are scaled to have a maximum circular 
velocity of approximately 60 km/s, which is the 
value for dwarf galaxies such as NGC3109 (Jobin 
& Carignan 1990) and NGC6822 (Weldrake, de 
Blok & Walter 2002). The virial masses for these 
models are in the range (4 — 7) x lO^^M© . In all 
cases the baryonic component represents a vciry 
small fraction (l%-2%) of the virial mass, which 
is a factor of ten smaller than the cosmological 
ratio of the baryons to the dark matter. 

All models have extended (r-^ 100 kpc) halos, 
which have high concentration NEW density pro- 
files. We did not try to make more realistic fits 
for the dark matter halos. In reality the halos 
should be slightly perturbed in the central part 
because of the presence of baryons. We decided 
to neglect this effect in favor of simplicity of the 
models. Most of our models have very thin disks. 
Real disks are typically thicker than what we as- 
sume for our models. The only model, which has 
a realistic disk is the "Disk" model V. 

The first two models (dubbed "Bar" and 
"Dwarf") start with a pure stellar disk embed- 
ded in an extended dark matter halo. Initially the 
disks in the models are submaximal: about 50% 
of mass in central 1 kpc is in the disk the rest is 
in the dark matter. The disks are relatively cold 
(stability parameter Q is 1.2). As the result of this 
combination of parameters, the disks are unstable 
and produce bars. The main difference between 
the models is that the disk in the model "Bar" 
is significantly more massive than in the model 
"Dwarf" . 

Models "Bulge", "Thin Disk", and "Disk" do 
not have bars. The models are stable - they do 
not evolve much. Nevertheless, we run the models 
for about 1.5-2 Gyrs to ensure detailed equilib- 
rium. Indeed, small variations are found in the 



7 



models. Initially the rotational velocity does not 
depend on the distance z to the disk plane. Af- 
ter evolution we find that the rotational decreases 
with z. This tendency increases effects of galaxy 
inclination. Models "Thin/Thick Disk" have rela- 
tively hot disks (Q = 3.0). The disks are also more 
extended, which make them sub-dominant even in 
the center. The disks did not develop a bar. 

Figure 4 presents details of the distribution of 
different quantities for models II and III. The mod- 
els to some degree are similar and are representa- 
tive for the other two models (I and IV). In mod- 
els II and III the dark matter is dominant in the 
outer regions, but is sub-dominant in the central 
500 pc. Velocity dispersions of baryonic compo- 
nents are very small. For example, for models II 
and III the rms velocities are smaller than 15 km/s 
even at Ikpc radius. 

5. The Numerical Simulations 

The procedure of generating of initial condi- 
tions is described in detail in Valenzuela, & Klypin 
(2003). We use the following approximation for 
the density of the stellar disk in cylindrical coor- 
dinates: 



pdiR,z) = 



e '^d sech^{z/zo), (1) 



where is the mass of the disk, Ra is the expo- 
nential length, and zq is the scale height. The 
later was assumed to be constant through the 
disk. Effective exponential scale height (height at 
which the density decreases by e— times) is equal 
to Ze = 1. 08520. The disk is truncated at five 
exponential lengths and at three scale heights: 
R < 5Rd, z < 3zo- Disk mass was slightly cor- 
rected to include the effects of truncation. The 
vertical velocity dispersion is related to the sur- 
face stellar density E and the disk height zq- 



(jI{R) = nGzoT.{R), 



(2) 



where G is the gravitational constant. The ra- 
dial velocity dispersion an is also assumed to be 
directly related to the surface density: 



4(i?) = Ae 



(3) 



The normalization constant A in eq. (3) is fixed 
in such a way that at any radius R the rms ra- 
dial random velocity is Q times the critical value 



needed to stabilize a differentially rotating disk 
against local perturbations. The tangential (</>) 
component of the rotational velocity and its dis- 
persion arc found using the asymmetric drift and 
the epicyclic approximations. For the disk compo- 
nent we use the thin-disk approximation. 

For model "Bulge" we also add a spherical non- 
rotating component with the Hernquist density 
profile: 



Pdm(r-) 



r/vH, 



(4) 



where th is the 1/4-mass radius. For the model 

Th = 125 pc, which gives half-mass radius 250 pc. 

We assume that the dark matter density profile 
is described by the NFW profile (Navarro, Frenk, 
& White 1997): 



/5dm ('f) 



Ps 



x{l+xY 



, X ^r/rs 



ln(l + C) - 



C 



l + C 



C = 



r. 



(5) 

(6) 
(7) 



where Myir and C are the virial mass and the 
concentration of the halo. For given virial mass 
the virial radius of the halo is found assuming a 
flat cosmological model with matter density pa- 
rameter SIq = 0.3 and the Hubble constant Hq = 
70km/sMpc"^ 

Knowing the mass distribution of the system 
M{r), we find the radial velocity dispersion of the 
dark matter: 



2 



1 

Pdn 



' GM{r) ^ 
Pdm 5 — dr.. 



(8) 



The other two components of the velocity disper- 
sion are equal to cr^dm- other words, the ve- 
locity distribution is isotropic, which is a good 
approximation for the central parts of the dark 
matter halos. At each radius we find the escape 
velocity and remove particles moving faster than 
the escape velocity. The dark matter halo is trun- 
cated at the virial radius, 

Velocities are picked from an appropriate Gaus- 
sian distribution truncated at the escape velocity. 
The disks are realized with particles of an equal 
mass. The dark matter halo is composed of parti- 
cles of variable mass with small particles placed in 
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the central region and larger particles at larger dis- 
tances. We double the particle mass as we move 
further and further from the center. The region 
covered by small mass particles is ^ 20 kpc much 
larger than the size of the disk. In order to reduce 
the two-body scattering each dark matter particle 
in the central region has the same mass as a disk 
particle. This procedure is designed to reduce the 
number of particles and, at the same time, to allow 
us to cover a very large volume. 

We use the Adaptive-Refinement-Tree (ART) 
A''-body code (Kravtsov et al. 1997; Kravtsov 
1999) to run the numerical simulations analyzed 
in this paper. The code starts with a uniform grid, 
which covers the whole computational box. This 
grid defines the lowest (zeroth) level of resolution 
of the simulation. The standard Particles- Mesh al- 
gorithm is used to compute the density and grav- 
itational potential on the zcroth-lcvcl mesh with 
periodical boundary conditions. The code then 
reaches high force resolution by refining all high 
density regions using an automated refinement al- 
gorithm. The refinements are recursive. A refined 
region can also be refined. Each subsequent re- 
finement level has half of the previous level's cell 
size. This creates a hierarchy of refinement meshes 
of different resolution, size, and geometry covering 
regions of interest. Because each individual cubic 
cell can be refined, the shape of the refinement 
mesh can be arbitrary and effectively match the 
geometry of the region of interest. This algorithm 
is well suited for simulations of a selected region 
within a large computational box, as in the simu- 
lations presented below. 

The criterion for refinement is the local density 
of particles. If the number of particles in a mesh 
cell (as estimated by the Cloud-In-Cell method) 
exceeds the level ntiucsh, the cell is split ("refined") 
into 8 cells of the next refinement level. The refine- 
ment threshold depends on the refinement level. 
The threshold for cell refinement was low on the 
zeroth level: nthrosh(O) = 2. Thus, every zeroth- 
levcl cell containing two or more particles was re- 
fined. The threshold was higher on deeper levels 
of refinement: rithresh = 3 and rithresh = 4 for the 
first level and higher levels, respectively. 

During the integration, spatial refinement is ac- 
companied by temporal refinement. Namely, each 
level of refinement is integrated with its own time 
step, which decreases by factor two with each re- 



finement. This variable time stepping is very im- 
portant for accuracy of the results. As the force 
resolution increases, more steps are needed to in- 
tegrate the trajectories accurately. 

6. Analysis 

We have analyzed the N-body simulations using 
several methods or algorithms that mimic those 
used by the observers. Our data consist of ve- 
locities, positions and masses of the particles that 
constitute the disk of the model galaxy. The sim- 
ulations are set up so that the plane of the disk 
is the X — y plane. During our 'observations' we 
'look' down the 2;— axis of the simulation. To ob- 
serve our galaxy at various inclinations we simply 
rotate y — z plane about the a;— axis by varying 
angles. 

The object of the exercise is to calculate rota- 
tion curves for our model galaxies as well as re- 
construct the density profile. We have also con- 
structed velocity field diagrams in the form of ve- 
locity contour plotted on the sky. One method 
frequently used by optical observers is long slit 
spectroscopy. In oxvc case, we place a rectangular 
slit running parallel to and centered on the x-axis. 
We then select all particles whose x—y coordinates 
place them within the slit boundaries. We define 
pixels along the slit and we assign particles falling 
in the slit to these pixels. We compute the line-of- 
sight velocity of each particle. We then compute 
the average velocity of the particles in each pixel. 
This average velocity, determined as a function of 
distance along the slit, is our "observed" rotation 
curve. By comparing this average rotation veloc- 
ity with the true average velocity as well as the 
circular velocity, we can estimate the accuracy of 
the method. 

In the case of HI observations in the radio and 
Fabry-Perot measurements in the optical, one has 
velocity information over the whole galaxy image 
not just along a narrow slit. In this case, observers 
determine the rotation curve using a tilted ring 
model of the velocity field. The method has been 
described by Begeman (1989). We reproduce the 
salient points below 

6.1. The Tilted Ring Model 

The model is based on the assumption that a 
rotating disk galaxy can be described by a set of 
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concentric rings. Each ring has a constant cir- 
cular velocity Vc and two orientation angles (the 
inclination angle, i, and the position angle of the 
major axis, (f>). The values of the ring parameters 
are determined from the observed radial velocities 
in a set of concentric elliptic annuli in the plane 
of the sky. The ellipticity of the annuli, deter- 
mined from isophotal analysis, fixes the inclina- 
tion of each ring. Vq is the systemic velocity of 
the galaxy (in our case km/s). For a given ring, 
the observed radial velocities recorded on a set of 
sky coordinates {x, y) are related to the basic pa- 
rameters by: 

V{x, y) = Vo + Vc{R) sin(i) cos{e). (9) 

In the expression above, 6 is the azimuthal angle in 
the plane of the galaxy, related to the parameters 
i,(f),Xo and yo through: 

^Qg^^^ ^ -(a: - Xq) sm{(j)) + {y - yp) cosjcj)) ^^^^ 

R 

dnfgl - ~ ^^^^"^^ ~ ^'"^^^ fill 

^ ' Rcos{i) ' ^ ' 

where R is the mean radius of the ring in the plane 
of the galaxy. 

The first step in such a procedure is to deter- 
mine the the position angle and ellipticity of the 
'isophotes' of the model galaxy image. To do this 
we project all the particle positions on the sky 
plane. We then form an image in FITS format 
that we process using IRAF ^. We use a package 
within IRAF called STSDAS 2. Within STSDAS 
we use the ELLIPSE task which to fit a set of el- 
liptical isophotes over the image. The method has 
been described by Jedrzejewski (1987). 

Line-of-sight velocities measured at positions 
close to the minor axis carry less information 
about the underlying circular velocity than those 
measured at positions close to the major axis. We 
have excluded all values of the angle 9 > 45° where 



^IRAF is written and supported by the IRAF programming 
group at the National Optical Astronomy Observatories 

(NOAO). 

^Tliis is the software for reducing and analyzing data from 
the Hubble Space Telescope. It is layered on top of IRAF 
and provides general purpose tools for astronomical data 
analysis as well as routines specifically designed for HST 
analysis. 



9 is the angle with respect to the major axis in the 
plane of the galaxy. Weldrake, de Blok & Wal- 
ter (2002) experimented with different weighting 
schemes, such as down weighting the minor axis by 
I cos^l and cos^ 9. They also varied the exclusion 
by angle between 45°and 15°and found no differ- 
ence in the curves produced for their analysis of 
the galaxy NGC 6822. 

Once we know the values of 4>, i and R we 
can compute Vc for the each elliptical anrnilus and 
thus obtain Vc{R), the rotation curve. As with 
the long slit method, we can compare Vc{R) with 
the known rotational velocities and true circular 
velocities for the model galaxy. 

6.2. Density profiles 

One can invert the observed rotation curve to 
determine the parent mass distribution (de Blok, 
McGaugh, Bosma, & Rubin 2001) using the equa- 
tion 



2vdv 

AnGpir) = - — 



(12) 



This procedure assumes that the galaxy is spher- 
ically symmetric, completely dark matter domi- 
nated and that the gas moves in circular orbits in 
a planar disk. It is worth checking the reliability 
of the method using a realistic model of a galaxy. 

7. Results 

7.1. Long Slit Observations 

7.1.1. Models IV and V (no bar, no bulge) 

We observe the models IV and V in a manner 
designed to mimic long slit observations. In the 
following discussion we take circular velocity to 
mean the velocity at which a particle would move 
if it was moving in uniform circular motion around 
the center of the galaxy. The term rotational ve- 
locity refers to the speed at which a particle actu- 
ally is rotating about the center of the galaxy; the 
total velocity of the particle can of course be larger 
than this. The term "observed velocity" refers to 
the mean observed velocity at a given radius from 
the center. As we shall show, this latter number 
can have a value that is considerably smaller than 
the rotational velocity which in turn can be con- 
siderably smaller than the circular velocity. 

Figure 5 compares the rotation curve obtained 
in this manner with true rotation velocities in the 
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disk which we can measure directly. We also de- 
termine the total circular velocity that would be 
measured if there were no random velocities in the 
disk (i.e. perfect circular motion in the galaxy 
potential). In order to mimic slit observations, 
we incline the model by some angle and find line- 
of-sight velocities for "stars" inside a long slit of 
100 pc width. 

A clear trend is that the discrepancy between 
the "observed" and true rotation curves increases 
with inclination of the galaxy. By the time one 
reaches an inclination angle of 80° the discrepancy 
is ~ 8 km/s at a radius of 0.5 kpc. This is a 30% 
in the estimation of velocity in the central part of 
the galaxy. 

Why does the long slit mode of observation lead 

to underestimates of the rotational velocities in 
the disk? Let us consider a simple example. For 
model IV we take a line of sight through the plane 
of the disk whose closest approach to the center 
of the galaxy is 0.7 kpc. We choose this number 
since it is within the central kiloparsec of dwarf 
galaxies that mismatch between theory and ob- 
servation is greatest. Averaging the velocities of 
N-body particles along this line of sight is equiva- 
lent to measiiring the rotation velocity at one pixel 
location of a long slit observation. The velocities 
are weighted by the density of material along the 
line of sight. Ironically, the smallest rotational ve- 
locity occurs for the points that arc located closest 
to the center of the galaxy. Since we are averag- 
ing this velocity with velocities of points further 
from the center that thus have larger circular ve- 
locities, how can the average be lower? The point 
of course is that we measure radial velocities along 
the line of sight. For example, a particle located 
5 kpc away from the center of our model rotates at 
60 km/s. However, the projected velocity of this 
particle on a line of sight 0.7 kpc from the center 
of the model is only 9 km/s. The rotational veloc- 
ity of the disk 0.7 kpc from the center is 30 km/s. 
The average velocity that we measure along this 
line of sight is 20.4 km/s or 70% of the true ve- 
locity. One might well expect this difference to 
be even larger given that we are including radial 
velocities as low as 9 km/s. The answer is that 
the integral is intensity weighted so that points 
near the center carry considerably more weight in 
the mean determination than those near the edge 
of the galaxy. For example, for the whole line of 



sight we get a mean velocity of 20.4 km/swhile we 
get 26 km/sif we include only the central kpc. 

The discussion above refers to the model viewed 
edge-on. The effect of viewing the galaxy at a 
lower inclination is to include particles that are 
above and below the galactic plane in the inte- 
gral. Figure 6 shows the radial velocity versus line 
of sight distance for particles viewed with an in- 
clination angle of 75°. The points in the plot on 
the left are selected at location along the slit of 
0.6 to 0.8 kpc, those in the right hand plot are 
selected from 4. to 4.2 kpc. It is clear (as ex- 
pected) that particles located above the plane of 
the galaxy have slower rotational velocities. The 
mean velocity measured as a function of distance 
along the line of sight matches the rotation curve 
computed for particles 0.25 kpc above the disk at 
about 1 kpc, which is where one would expect it 
to for an inclination angle of 75°. Note that par- 
ticles further than 0.25 kpc from the plane of the 
galaxy are not substantially affecting the velocity 
measurement because the density out there is too 
low. 

One might ask why, given these projection ef- 
fects, the rotational velocity is only underesti- 
mated for lines of sight passing close to the center 
of the galaxy. Lines of sight that pass more than 
-^2 kpc from the center the measured velocities 
are in good agreement with the true rotational ve- 
locities for our models. The right hand plot of 
Figure 6 illustrates this; we do a good job recover- 
ing the rotational velocity of the particles at large 
radii. When the distance of the line of sight from 
the center of the galaxy is comparable to the size 
of the disk, the projection effects are quite differ- 
ent to those discussed in the previous case. Firstly, 
the spread in radii of particles in the line of sight 
is smaller at for lines of sight at larger distances 
from the center. Secondly, the projection effects 
are less pronounced at larger radii. Particles going 
at a given rotational speed have higher measured 
radial velocities because of the inherent geometry. 
These points are illustrated in Figure 7. 

1.1.2. Model I: The Barred Galaxy Model 

We refer to Model I in Table I as the barred 
galaxy model. For the case of this model, the sit- 
uation is much more complex. We have examined 
3 cases which span the range of possible geome- 
tries. One with the bar at 45 degrees to the line of 
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sight (referred to as bar45), one case with the bar 
ahgned along the hne of sight (barO) and thirdly a 
case where the bar is aligned perpendicular to the 
line of sight (bar90) . In each case wc use the same 
galaxy model; we just change the viewing angle. 
For each of these three cases we choose different 
inclination angles. 

In the bar45 geometry viewed at an inclination 
of 60° we recover the true rotation curve with an 
accmacy 2 km/s. In the case of a bar aligned with 
the line of sight (barO) one gets larger rotational 
velocities measured in the central kpc region of 
the galaxy than are truly present. This is due 
to motions within the bar that are not rotational 
but contribute to the perceived rotational velocity. 
Such a bar viewed lengthwise would not be classi- 
fied as a bar but most likely as bulge and likely will 
be discarded. The velocity differences observed in 
this situation can be large; ~ 10 — 15km/s. The 
third case when the bar is perpendicular to the 
line of sight and clearly visible as a bar, results 
in observed rotational velocities that are smaller 
than the true rotational velocities in the disk. The 
mismatch between the long slit result and the ro- 
tational velocities for a galaxy viewed at an incli- 

o 

nation angle of 75 is of about 25% at a 1 kpc 
radius. Inclination angle determinations are sub- 
ject to rather large uncertainties as we mentioned 
above. If an inclination angle of 65° had been 
determined this galaxy the mismatch between the 
true and observed velocities would be increased by 
several km/s. As we showed in Figure 3 it is galax- 
ies with bars that tend to have the lower rotational 
velocities in their centers. 

The conchision from the long slit analysis is 
that it is possible to underestimate the rotational 
motions in the central parts of dwarf galaxies by 
as much as 10 km/s using long slit data. This is 
a 30% error in the rotational velocity. The error 
is simply due to the way the line of sight samples 
the velocities in a rotating disk of finite thickness. 
Adding a bar into the picture can enhance or de- 
crease this effect depending on the orientation of 
the bar relative to the observer. When the bar 
lies in the plane of the sky the differences are in- 
creased. Bars are not always easy to detect. How- 
ever the inclination angle and position angle mea- 
sured as a function of radius can give indications 
as to the presence of a bar. 



7.2. The Tilted Ring Model: Results 

1.2.1. Model IV (no bar, no bulge) 

At low inclinations we find that a small error in 
the inclination measurement can affect the veloc- 
ity results substantially. For example, we viewed 
model IV with at an inclination of 10° but the 
ellipse fitting software estimated an inclination 
14° resulting in underestimates of the rotational 
velocity of 10-15 km/s about 40%. At inclinations 
of ~ 60° the measured velocities accurately match 
the rotational velocities over the measured range 
of 0.2 to 4 kpc. At an inclination of 75°wc find that 
underestimates of the velocity of ~ 5km/s occur 
at radii of less than 0.4 kpc, a 20% error in veloc- 
ity. At an inclination of 80° the velocity under- 
estimates occur at radii smaller than 0.7 kpc. At 
0.2 kpc the mismatch is about 8km/s, a 25% error 
in the velocity estimate. These results are shown 
in figure 8. 

7.2.2. Model I: Effect of a bar on the tilted ring 

results 

The tilted ring analysis described above is the 
one adopted to calculate rotation curves using well 
resolved high signal to noise data. The method 
is applied to datasets where velocity information 
is available over the whole image plane, not just 
along a narrow slit. In practice observers have 
assumed a fixed or slowly varying position angle 
and inclination as a function of radius when imple- 
menting this method (Weldrake, de Blok & Walter 
2002). 

To see the effect of the bar on rotation velocity, 
we selected disk particles in a ring from 1 kpc to 
1.1 kpc and examined their rotation velocity as a 
function of azimuthal angle. The signature of the 
bar is clear. Particles located close to the position 
angle of the bar have significantly smaller rota- 
tional velocities on average. The bar is located 
two quadrants of the galaxy. We compute the av- 
erage velocity of all the particles in the ring and 
obtained a value of 38.6km/s. When we exclude 
the quadrants containing the bar we obtain a value 
of 43.5 km/s. We have plotted the mean rotation 
velocities of particles in the disk, computed by ex- 
cluding particles in the bar, this is the dotted line 
shown in Figure 9. We also show as a dashed line 
the circular velocity in Figure 9. The point is that 
there are a number of significant mismatches in the 



12 



center of the galaxy. The tihcd ring method does 
not recover the average rotational velocities in the 
disk. Even if it did, this would not be sufficient, 
since the average is lowered by the presence of a 
bar. Even if one has a good measure of the rota- 
tional velocities in the disk, one needs an accurate 
method of calculating asymmetric drift to recover 
the circular velocity. The bottom line is that large 
mismatches occur between the observed and true 
velocities in the central part of the disk when a bar 
is present. At 1 kpc from the center of the galaxy 
the observed velocity is 30% less than the true ve- 
locity. At 0.5 kpc from the center of the galaxy 
the observed velocity is 50% of the true velocity. 

It does not take a very large error in velocity 
to produce a significant error in the power law 
slope of p{r). Let us suppose we compute p{r) 
at 1 kpc and 2 kpc from the center of the galaxy 
and that wo got the correct answer at 2 kpc but 
underestimate the rotation velocity and 1 kpc. Let 
us say we think the velocity is 31km/s when the 
true velocity is 35kni/s. The percentage error is ^ 
10%. The rotational velocity enters as the square 
in the mass measurement so the resulting error in 
the density is ~ 20% such that a power law index 
-1 is measured to be -0.65. 

7.2.3. Model I: Recognizing a Bar 

CCD imaging studies reveal that '--^70% of spi- 
ral galaxies are barred (Eskridge et al. 2002). To 
avoid complications, galaxies should be selected 
for rotation curve studies from the remaining 30%. 
However, our models suggest that bars may well be 
present in galaxies where imaging does not reveal 
their presence. We have demonstrated that bars 
in galactic disks give rise to errors when inferring 
galaxy masses from rotation curves. One might 
argue that barred galaxies are easily recognizable 
and hence can be easily excluded from galaxy sam- 
ples. This is not the case. We show in figures 10 
and 11 contour plots of the surface density of par- 
ticles for the barred galaxy model viewed at an in- 
clination of 75°. The surface density is measured 
as magnitudes and the contours are equally spaced 
in magnitude to mimic observational plots. It is 
certainly difficult to see the presence of a bar from 
density contours alone. Bars do cause variations 
in the inclination angle and position angles mea- 
sured as a function of radius but warps and bulges 
can produce similar effects. We also show next to 



the contour plots, the "luminosity profiles" of our 
model disks. These plots also reveal the presence 
of the bar in an indirect way. 

It is striking to what extent the signature of 
the bar can vary in the surface density profile de- 
pending on the viewing angle. The bar can have 
the appearance of a bulge, or produce a sudden 
decrease in surface brightness, or have no marked 
effect. 

7.2.4- Model ILEvolution of a bar 

We use model II to study effects of bar evolution 
on rotation curves. We have examined outputs of 
the simulation at three stages of the bar forma- 
tion. The first snapshot was made 1.4 Gyr after 
the start of the simulation when the bar is still 
weak. The second snapshot was made 1.62 Gyr 
after the start of the simulation at which point 
the bar is strong. The final snapshot was taken 
2.15 Gyr after the start of the simulation when 
the bar has buckled. Buckling refers to the bend- 
ing of the bar into a peanut shape forming a thick- 
ened double lobcd structure (Combes & Sanders 
1981). The rotation curve results are shown in 
Figure 12. As expected there is a large increase in 
mismatch between the circular and observed ve- 
locities as the strength of the bar increases. After 
the bar has buckled, the observed velocity is less 
than 50% of the circular velocity 0.5 kpc from the 
center of the galaxy. Figure 13 shows the incli- 
nation based on the isophote analysis. One can 
clearly see the effect of the bar extending to larger 
radii as the simulation evolves. Since we are view- 
ing the bar perpendicular to the line of sight there 
are no isophote twists although these would be 
present if one viewed the bar at some angle either 
than perpendicular or parallel to the line of sight. 

7.2.5. Model III: Effect of a bulge on the tilted 
ring results 

Infrared studies of spiral galaxies and low sur- 
face brightness galaxies (Boker, Stanek & van der 
Marel 2003; Galaz, Dalcanton, Infante, & Treis- 
ter 2002; Beijersbergen, de Blok & van der Hulst 
1999) show that the majority of these systems con- 
tain bulges. We have run a simulation to explore 
what effect the presence of a bulge has on the 
tilted ring results. Model III "Bulge" is designed 
to mimic the bulges detected by Boker, Stanek & 
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van dor Marcl (2003): it has a very small bulge 
with half mass radius is 250 pc. Figure 14 shows 
the results of the tilted ring analysis for model III. 
It is clear that bulges as well as bars can result 
in low rotational measurements in the central re- 
gions of galaxies where the circular velocity is in 
fact high. The differences are of order 10-20km/s 
in the central 0.5 kpc. The velocity error ranges 
from ~20% at 1 kpc to ~ 50% at 0.2 kpc. In the 
case of a barred model the greatest cause of the 
mismatch are the non-circular motions produced 
by the bar. Due to the presence of the bar, a large 
fraction of stars and gas parcels have velocities 
orthogonal to the line of sight. These velocities 
are undetected by the emission line observations. 
This effect strongly depends on the bar orienta- 
tion and is maximum when the bar is perpendic- 
ular to the line of sight. The case of galaxies with 
small bulges is slightly different. A bulge is a sys- 
tem substantially supported by random motions. 
An observer performed an average along the line 
of sight taking into account all particles belong- 
ing to either the disk or the bulge. For disk par- 
ticles the detected radial velocity is close to the 
projected circular velocity along the line of sight, 
a small correction should be applied because of 
the asymmetric drift. For bulge particles the de- 
tected radial velocity component is comparable to 
the velocity components orthogonal to the line of 
sight. Systematically bulge particles show smaller 
velocity along the line of sight compared with disk 
particles, for the same circular velocity. There- 
fore the avenge or the observed rotation velocity 
is artificially biased toward lower values. 

When models are viewed at low inclinations, 
error in the inclinations can cause large system- 
atic errors in the rotation velocity estimates. For 
example the tilted ring algorithm produces an in- 
clination value of 26.7° for the model viewed at 
an inclination of 20°. The resulting error in sini 
results in a ~25% error in velocity. A 25% error 
in velocity translates as a 50% error in the density 
in the outer parts of the galaxy. 

The dynamical signatures of the bulge model 
extend to larger radii than its signature in the 
surface density profile plots. For an inclination 
of 60°, the small bulge does not make a significant 
contribution to the surface density profile beyond 
400pc, but rotation velocities are determined to 
be significantly lower due to its presence beyond 



1.2 kpc. By comparison, the model with no bar 
and no bulge, model IV, when viewed at an incli- 
nation of 60° shows good agreement between the 
rotation velocity and the observed velocity down 
to 200 pc. 

7.2.6. Model V: Effect of disk thickness on the 

tilted ring results 

The motivation for Model V "Disk" is to gauge 
the effect of disk thickness on the underestima- 
tion of rotational velocity in the central Kpc of 
the models. Figure 15 shows velocity and inclina- 
tion for the disk model viewed at an inclination 
of 75°. The mismatch between the "observed" ve- 
locity and the rotational velocity is 10 km/s at 
0.4 Kpc and 6 km/s at 0.2 Kpc. The thin disk 
model IV analysis reveals a smaller mismatch of 
3 km/s at 0.4 Kpc and 6 km/s at 0.2 Kpc. As 
expected the mismatch increases with disk thick- 
ness at the largest inclinations. 

7.2.7. Surface Density Profiles 

Using the Ellipse task in the package STSDAS 

we can produce surface density profiles for our 
models. This is the same package used by Boker, 
Stanek & van der Marel (2003) among others to 
produce surface density profiles of galaxies. We 
are thus processing our A^— body "data" in as simi- 
lar a manner as possible to the observers. The sur- 
face density profile for the barred galaxy model is 
show in the four panels of Figure 16. We are view- 
ing the galaxy at an inclination of 75° and showing 
how the surface density profile is affected by vary- 
ing the angle the that bar makes with the line of 
sight. The bottom right panel showing the profile 
that results in viewing the bar along the line of 
sight, illustrates the point that in this geometry a 
bar can be mistaken for a bulge. When the bar 
is viewed perpendicular to the line of sight we see 
a sudden dip in the surface brightness at about 
2 kpc. This dip is not so prominent when the 
bar is only slightly rotated from the perpendicu- 
lar (upper right panel). It is clear that a number 
of the features found in these surface brightness 
plots qualitatively resemble features found in the 
data presented by Boker, Stanek & van der Marel 
(2003). 

Figure 17 shows the same analysis for model III 
(bulge model). By varying the inclination angle 
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wc vary the contribution of the disk to the surface 
brightness. The slope of the disk contribution also 
varies depending on the viewing angle. The bulge 
signature in the surface brightness plot docs not 
extend beyond ~ 0.6 kpc whereas the rotation ve- 
locity is underestimated out to ~ 1.2 kpc. 

Figure 18 shows the effect of viewing the model 
at constant inclination during its evolution. One 
can see the growth of the bar as a shift in the loca- 
tion of the point at whic;h the; profile slop changes. 
This point occurs at roughly 0.5 kpc in the upper 
panel, shifts to 1.5 kpc in the lower left panel and 
back inward to about 1 kpc as the bar buckles. 

Figures 16 to 18 illustrate the variety of profiles 
that can be produced by adding bulges and bars 
to rotating disks. In few cases, based on surface 
density profiles alone, is it immediately clear that 
a bar or bulge must be present. Bars can appear 
as bulges or disappear almost entirely depending 
on the viewing angle. Bulges can affect only the 
central 200 pc or extend as far as 1 kpc depend- 
ing on the intrinsic bulge size and viewing angle. 
The bar signature also changes as the instability 
evolves in the disk. We have shown earlier that 
bars and bulges significantly affect the interpre- 
tation of rotation curves. The figures discussed 
in this section show that is is difficult to unam- 
biguously identify these features based on surface 
density profiles alone. This section is intended as 
a cautionary note to researchers who ignore these 
features in the analysis of galaxy rotation curves. 

7.3. Density Profiles From Rotation Curves 

As explained in section 5.2, one can infer the 
density distribution in a galaxy from its rotation 
curve using eq. 12. Under the assumption that 
gas and stellar mass contribution is negligible, it 
is possible to recover the dark matter halo density 
profile. LSB and Dwarf galaxies are considered to 
closely satisfy the assumption of being completely 
dark matter dominated and it is the application 
of this is the method to LSB and Dwarf galaxies 
which gives rise to the inference of flat cores in the 
central density distributions of galaxies. With our 
models we can check the accuracy of this method. 

For model IV (with no bar and no bulge) the 
ellipse fitting method does a reasonable job recov- 
ering the true density. The density versus radius 
is shown for the model with no bar viewed at an 



inclination of 80° in Figure 19. At radii of less 
than 1 kpc the density is underestimated by about 
30%. The underestimates vary from 10% to 50% 
depending on the inclination angle. 

These effects are much more pronounced for the 
barred galaxy model. In figure 20 we show the 
density for the barred model viewed at an inclina- 
tion of 75°. Figure 21 and Figure 22 show the same 
analysis with the bar in the plane of the sky and 
perpendicular to the plane of the sky. The effects 
are most pronounced on the latter case. When the 
bar is perpendicular to the line of sight it is easy 
to see how an observer would infer the presence of 
a constant density core in the center of the galaxy. 
The orientation of the bar with respect to the line 
of sight has a strong effect on the inferred density 
profile. 

We have computed the the slope of the p{r) 
curve at locations smaller than 1 kpc for the three 
barred models. We show these in the histogram in 
Figure 23. The key points is that for barred galaxy 
models viewed at an inclination greater than or 
equal to 60°we always measure a slope that is shal- 
lower than the true slope. For a true slope of -2 
we can obtain slopes as shallow as 0. Note that 
the NFW profile has a slope of -1, in the light of 
our findings it is not surprising that slopes of shal- 
lower than -1 are being found near the centers of 
dwarf galaxies. 

8. Conclusion 

The rotation curves of dwarf and LSB galax- 
ies arguably present the most challenging prob- 
lem for current cosmological models. The key to 
this problem lies in reliably estimating the accu- 
racy with which one can recover the mass distri- 
bution and, thus, the dark matter density using 
observed rotation curves. This is the main ques- 
tion which we address in this paper by studying 
realistic numerical models. The models mimic real 
galaxies in many respects. For example, they have 
thin 100 pc) and cold 10 km/s) exponen- 
tial disks. The latter condition is very important 
because many effects are related to the magnitude 
of random motions. 

The main lesson to be learned from our analy- 
sis is that there are many systematic biases. Each 
bias can be very small (only few km/s), yet they 
all add up to produce significant effects. We find 
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that the tilted ring analysis of spectroseopie data 
results in large 30% to 50% underestimates of the 
rotation velocity in the central 1 kpc region of the 
galaxy models. Those underestimates in velocity 
in turn result in the inference of slopes in the cen- 
tral density distributions of our models that are 
shallower than the true slopes. 

A quantitative understanding of these biases is 
important if one is to accurately interpret rotation 
curves. If one ignores these biases, the analysis 
of rotation curves can produce misleading results. 
The extent to which this has occurred for real 
galaxies can only be resolved by detailed analyses 
of individual objects. One cannot be sure that the 
analysis of real observations data was done prop- 
erly when the same analysis produces false results 
when applied to realistic models. 

One of the difficulties in interpreting observa- 
tional data lies in the fact that galaxies are treated 
in an overly simplified manner (e.g., cold gas inside 
thin exponential disks). Observational results (in- 
cluding those presented in this paper) suggest that 
this is far too naive a picture. In the central region, 
the stellar component often rotates as fast as the 
emission-line gas - in contradiction with the pre- 
diction of a naive picture of a cold gas embedded 
in a hot stellar component. The presence small 
(200-500 pc) bulges and bars, which can be easily 
overlooked, further complicates the interpretation 
of the data. 
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Fig. 1. — Stellar and gas rotation velocities of 
observed galaxies. Symbols show velocities mea- 
sured using emission lines. The full curves with 
error bars are for absorption lines. Galaxies on 
the left panels have similar rotation in gas and in 
stars. For galaxies on the right gas rotates typi- 
cally faster than stars. Of those three galaxies are 
barred. 
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Fig. 2. — Surface brightness profiles. Dashed 
curves show disk/bulge double exponential ap- 
proximations. In the case of barred galaxies the 
exponentials do not join though the outer disk are 
still clearly exponential. The spike at 10 kpc in 
profile of NGC5746 is due to a star. The ring of 
NGC5389 is seen at 6kpc as a hump overlaid on 
exponential disk. 
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Fig. 3. — Inclination angle versus inner power law 
slope a for the galaxies in the de Blok, McGaugh, 
Bosma, & Rubin (2001) sample. The circled 
points correspond to galaxies that have been iden- 
tified by McGaugh, Rubin, & de Blok (2001) as 
being barred galaxies. There should be no correla- 
tion of slopes with the inclination. Lack of slopes 
at left top corner clearly indicates significant bi- 
ases in the sample. 




R(Kpc) R(Kp<!) 



Fig. 4. — Profiles for the final moments of the 
model II ( "Dwarf' , left panels) and the model III 
("Bulge", right panels). Lower panels show circu- 
lar velocities of different components in the mod- 
els. The dotted curves show the rotation velocity 
of the disk components. For both models the bary- 
onic component (disk+bulge) dominates inside the 
central 700 pc; radius. Middle panels show velocity 
dispersions for the baryonic components. The dot- 
dashed curves are for the vertical (z) velocity. The 
full and dashed curves are for tangential (in the 
plain of the disk) and radial velocities correspond- 
ingly. Note that velocities are very small. Only 
inside central 500 pc do they go above 20 km/s. 
The top panels present the surface density profiles 
of the baryonic components scaled to magnitude 
units with arbitrary zero points. For the barred 
"Dwarf" the full curve is for the surface brightness 
along the bar and the dashed curve is perpendic- 
ular to the bar. The length of the bar is 1.8 kpc. 
The full curve for the "Bulge" model shows the 
total disk-|-bulge surface brightness. Disk compo- 
nent is shown by the dashed curve. 
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Fig. 5. — Velocities in Model V. The bottom 
panel shows the true rotation velocity (full curve) 
and the circular velocity (dotted curve) for this 
model. Long slit observations underestimate the 
rotational and circular velocities. The effect in- 
creases with inclination angle. The dot-dashed 
(long and short dashed) curves are for inclination 
70° (80° and 90°). The top panel presents rms 
velocities of the stellar component. The random 
velocities are small (10 -15km/s), but they are 
responsible for the difference between the (true) 
rotation and the circular velocities. 



Fig. 6. — Distribution of line-of-sight velocities 
versus distance along the line of sight for particles 
selected at 0.6 to 0.8 kpc (left panel) and at 4.0 to 
4.2 kpc (right panel) along a " slit", centered on 
the major axis for Model IV "Thin Disk" viewed 
at an inclination of 75°. The solid line shows the 
average velocity of particles within 0.25 kpc of the 
galactic plane, projected along our chosen line of 
sight. The dashed line shows the projected rota- 
tional velocities of particles further than 0.25 kpc 
from the galactic plane. The points with error bars 
show the mean velocity in bins 0.1 kpc wide. 
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Fig. 7. — Radial velocities along a "slit" for model 
IV inclined by 75°. All velocities are corrected for 
the inclination. Bottom panel: distribution of ve- 
locities of individual particles. The curve shows 
the true rotation velocity. The top panel shows 
the ratio of observed to circular velocities. The ve- 
locity is underestimated by 20% at 0.5 kpc, while 
at distances larger than 2 kpc the errors are very 
small. 
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Fig. 8. — Rotation curves for model IV, the model 
with no bar and no bulge. The solid line shows 
the true circular velocity. The dashed line shows 
the rotational velocity. The points show the "ob- 
served" velocities for the model viewed at inclina- 
tions of 60°, 75°and 80°. The effects of inclina- 
tion are very small for 60°and for radii larger than 
1 kpc. The mismatch between the observed veloc- 
ities and the true velocities at small radii increases 
with increasing inclination. 
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Fig. 9. — Solid line in the bottom panel shows 
the average rotation velocity of particles in the 
disk as a function of radius. The dotted line is the 
same as the solid line but particles belonging to the 
bar have been excluded from the calculation. The 
dashed line is the circular velocity of particles in 
the disk. This is the velocity that particles would 
have if they moved in pure circular motion in the 
galaxy potential. At 0.5 kpc; the measured velocity 
is 40% of the true circular velocity. 
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Fig. 10. — Left hand side shows contour plot of the 
surface density of particles in the barred galaxy 
viewed at an inclination of 75°. The bar is per- 
pendicular to the line of sight. Right hand side 
shows surface brightness versus semi-major axis 
length for the same model. The surface brightness 
is measured on a magnitude scale. Presence of the 
bar is hardly noticeable in the contours of bright- 
ness. In the surface brightness the bars shows as 
a 0.5 mag feature at 2 kpc. 
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Fig. 11. — Same as Figure 10, except the az- 

imuthal angle of the bar is 45°resulting in twist- 
ing isophotes. The bar is not visible in the sur- 
face brightness profile. Presence of the bar is indi- 
cated only by the change in the position angle of 
isophots. 
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Fig. 12. — Evolution with time of the mismatch 
between the circular velocity and the observed 
rotation velocity. The mismatch gets larger as 
bar grows. Circular velocity (solid line), rota- 
tional velocity (dashed line) and 'observed veloc- 
ity' (points) for the dwarf galaxy model viewed 
during the formation and buckling of the bar. The 
model is viewed at an inclination of 60°. The bar 
is oriented perpendicular to the line of sight. 
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Fig. 13. — Inclination angle as a function of ra- 
dius for the dwarf galaxy model viewed during the 
formation and buckling of the bar. The model is 
viewed at an inclination of 60°. The bar is ori- 
ented perpendicular to the line of sight. The solid 
line shows the true inclination of the galaxy. The 
amplitude of the variation in the inclination an- 
gle does not change with time, but the radius of 
the transition to the true value increases with time 
and is a good measure of the bar length. 
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radius (Kpc) i = 60° 

Fig. 14. — Rotational velocity (bottom panel) and 
inclination angle (top panel) as a function of radius 
for model III (bulge) viewed at an inclination angle 
of 60°. The solid line in the bottom panel shows 
the circular velocity of particles in the disk. The 
dashed line shows the rotational velocity and the 
points show the observed velocity. The solid line 
in the upper panel is the mean inclination angle 
computed as the mean evaluated between 4 and 
5 kpc. 
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Fig. 15. — Model V "Disk" observed at an inclina- 
tion of 75°. The solid line shows the true circular 
velocity. The dashed line shows the rotational ve- 
locity. The points show the "observed" velocities 
for the model. The effects of inclination are very 
small for 60°and for radii larger than 1 kpc. The 
mismatch between the observed velocities and the 
true velocities at small radii increases with increas- 
ing inclination. 
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Fig. 16. — The four panels show the barred galaxy 
model viewed at an inclination of 75°for varying 
angle that the bar makes with respect to the line 
of sight. In the upper left panel the bar is viewed 
perpendicular to the line of sight. In the bottom 
right panel the bar is viewed along the line of sight. 
The two other panels are intermediate cases. 
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17. — Surface brightness profiles for model 
III (bulge model) viewed at an inclination of 
20° (dotted line) 60° (long dashed line) 80° (short 
dashed line). The solid line shows the surface 
brightness profile of the bulge on its own. The 
bulge contribution to the observed profile gets 
smaller as the inclination increases. Surprisingly, 
such a small bulge affects the rotation curve and 
the measured inclination angle up to 1 kpc as in- 
dicated in figure 14. 



25 



J i 1 i 1 1 i 1 i 1 1 1 1 1 1 1 1 1 1. 
^ 1.4Gyr 

7 \ -. 




: 1 1 1 1 1 1 1 1 1 1 M 1 1 1 1 1 i: 
^ 1.62Gyr - 

i « - 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1. 

K 2.15Gyr - 



+ i I I I I I 

12 3 12 3 



Radius kpc 



Fig. 18. — Dwarf galaxy model viewed at varying 
times during its evolution. 1.4 Gyr corresponds 
to the beginning of bar formation. By 1.62 Gyr 
a strong bar has formed, at 2.15 Gyr the bar has 
buckled as described in the text. The bar is viewed 
at 45°to the line of sight in each instance 
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Fig. 19. — Density versus radius for model IV (no 
bar, no bulge) viewed at an inclination angle of 
80°. The solid line is the actual density of the 
model. The points show the observed density and 
the dashed line is a fit to the observed density. 



Fig. 20. — Density versus radius for model I 
(barred model) viewed at an inclination angle of 
75°. The solid line is the actual density of the 
model. The points show the observed density and 
the dashed line is a fit to the observed density. The 
orientation of the bar is 45°to the line of sight. 




Fig. 21. — Density versus radius for model I 
(barred model) viewed at an inclination angle of 
75°. The solid line is the actual density of the 
model. The points show the observed density and 
the dashed line is a fit to the observed density. 
The bar is perpendicular to the line of sight. 
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Fig. 22. — Density versus radius for the barred 
model viewed at an inclination angle of 75°. The 
solid line is the actual density of the model. The 
points show the observed density and the dashed 
line is a fit to the observed density. The bar is 
viewed along the line of sight. 
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Fig. 23. — Histogram of the slopes obtained from 
model I (bar) viewed at inclination angles of 80 75 
and 60 °. The slopes where measured at positions 
evenly spaced between 0.2 and 1 kpc. The true 
slope of the model is -2. 



Table 1 

Parameters of galaxies with measured rotation curves 



Galaxy 


Mb 


Vlg 


Type 


Incl. 


-Rdisk/ -Rbulge 


Scale 


Comments 






km/s 




deg 


kpc 


kpc/" 




IC 267 


-19.91 


3577 


SBb 


45.0 


2.0/ 0.4 


0.248 


Very long 10 kpc bar. Compact 
















bulge. 


IC 1029 


-21.00 


2520 


Sb 


76.0 


5.0/ 0.4 


0.175 


Compact bulge. Gas rotation is 
















mcasurcid using SII emission line. 


NGC 5389 


-19.17 


1996 


SABa 


78.5 


3.0/ 0.45 


0.14 


Strong ring at 6 kpc. 


NGC 5746 


-20.51 


1676 


SABb 


80.5 


2.5/ 1.2 


0.12 


Large peanut-shape bulge/bar 


NGC 5719 


-18.79 


1676 


SABab 


70.0 


2.2/ 1.1 


0.116 


Barred galaxy with a large bulge. 


NGC 5793 


-19.11 


3387 


Sb 


70.5 


4.2/ 1.1 


0.234 




NGC 6368 


-19.98 


2904 


Sb 


75.5 


6.0/ 0.7 


0.201 


Gas rotation is measured using 
















SII emission line. 


NGC 7537 


-9.42 


2888 


Sbc 


70.0 


5.5/ 1.4 


0.200 
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Table 2 
Parameters of the Models 



Parameter Model 


Model I 
"Bar" 


Model II 
"Dwarf" 


Model III 
"Bulge" 


Model IV 
"Thin Disk" 


Model V 
"Disk" 


Halo mass Mhaio,MQ 


4.3x10^" 


6.0x10^" 


6.8x10^" 


6.0x10^" 


6.0x10^" 


Disk mass MdiskiM© 


9.5x10* 


5.7x10* 


6.8x10* 


5.7x10* 


5.7x10* 


Bulge mass Mbuige, Mq 






8.2 xlO^ 






Initial disk scale length R^, kpc 


1.16 


0.7 


1.75 


2.14 


1.75 


Final disk scale length R4, kpc 


1.50 


0.7 


1.75 


2.14 


1.75 


Initial disk scale height Zd, pc 


90 


108 


55 


150 


300 


Final disk scale height Zd, pc 


135 


120 


70 


150 


300 


Disk stability parameter Q 


1.2 


1.2 


1.5 


3.0 


3.0 


Bulge half-mass radius, pc 






250 






Number of halo particles 


3.3x10*^ 


3.45x10*5 


2.2 xlO^ 


3.45x10^ 


3.45 xlO^ 


Number of disk particles 


2 xlO^ 


1.6x10^ 


8x10^ 


1.6x10^ 


1.6x10^ 


Number of bulge particles 






9.7 xlO^ 






Halo concentration C 


11.5 


12.0 


11.0 


12.0 


12.0 


Maximum Force Resolution, pc 


11 


16 


19 


16 


16 


Duration of evolution, Gyrs 


4.2 


2.2 


1.4 


1.6 


1.2 


Time-step, yrs 


4 X 10^ 


1.4 x 10^ 


9.1 x 10^ 


1.4 X 10^ 


1.4 x 10* 
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